Rapidity Dependence of J/tp Production at RHIC and LHC 
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The motion of charmonium in heavy ion collisions is described by a three dimensional transport 
equation with initial production and continuous regeneration in hot medium. The observation 
of apparently stronger J/tp suppression at forward rapidity compared to that at midrapidity, so 
called J/tp puzzle at RHIC, can well be explained by the competition between the two production 
mechanisms. At LHC, however, the rapidity dependence of the J/tp production is dominated by the 
regeneration process. 
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PACS numbers: 25.75.-q, 12.38. Mh, 24.85.+p 



From lattice quantum chromodynamics (QCD) calcu- 
lations, a new state of matter, the so-called quark-gluon 
plasma (QGP), has been predicted to exist at high tem- 
perature and/or high baryon density. The J/tp produc- 
tion has long been considered as a probe of the QGP 
formation in relativistic heavy ion collisions lj. Differ- 
ent from the J / tp suppression observed at the CERN Su- 
per Proton Synchrotron (SPS) where almost all the char- 
monia arc primordially produced through hard nucleon- 
nucleon collisions (NN) and then suffer from nuclear ab- 
sorption 0, Q in the initial stage and anomalous sup- 
pression in the hot and dense medium 0, [f| @, 0, H, [§] in 
the later stage, there are a remarkable number of charm 
quarks in the QGP produced in high energy nuclear colli- 
sions at the BNL Relativistic Heavy Ion Collider (RHIC) 
and the CERN Large Hadron Collider (LHC), and the 
regeneration, namely the continuous recombination of 
those uncorrelated charm quarks offers another source 
for J/ip production 0, QH, QJ, d d • Obviously, the 
regeneration will enhance the J/tp yield and alter its mo- 
mentum spectrum. At midrapidity, the initial produc- 
tion and regeneration are almost equally important and 
the competition between the two production mechanisms 
controls [HI, [l6| the J/tp nuclear modification factor Raa 
and averaged transverse momentum square (p 2 ) observed 
at RHIC UlEl. 



Recently, the rapidity dependence of J/tp production 
in Au+Au collisions is measured at RHIC [lj], [3 and 
discussed in models [2(| Hi], [22j. As one can see, as a 
function of the number of participant nucleons N p , the 
observed Raa{Np) shows that the suppression is signifi- 
cantly stronger at forward rapidity (|y| <G [1.2,2.2]) than 
that at midrapidity (\y\ < 0.35). The enhanced suppres- 
sion is especially clear in central collisions, see Figs. 4 and 
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is about 0.6. In addition, we observe that the value of 
(p 2 ) is lower at midrapidity than that at forward rapid- 
ity, see Table II of Ref. [17j. The above observation is 
difficult to be explained by models with only initial char- 
monium production 17[ ■ Since the medium temperature 



at midrapidity is higher than or at least the same as that 
at forward rapidity, the suppression at forward rapidity 
is predicted in these models to be less than or at most 



the same with that at midrapidity. 

In this Letter, we investigate the rapidity dependence 
of J/tp production by constructing a three dimensional 
transport equation for the charmonium motion in the 
hot medium, with both initial production and regenera- 
tion. When the regeneration mechanism is included, the 
above J/tp puzzle is possible to be explained [2l|. Since 
the charm quarks are mainly distributed in the central 
region, their recombination into J/tps at midrapidity is 
more important than that at forward rapidity. While the 
medium absorption at midrapidity is in principle larger 
than that at forward rapidity, the competition between 
the regeneration and absorption may lead to an enhance- 
ment in Raa at midrapidity. We note that the effect in 
the rapidity dependence of Raa should also be elabo- 
rated with the rapidity dependence of (p 2 ). Considering 
the fact that the measured forward rapidity is still located 
at the plateau of light hadron rapidity distribution [23| , 
the temperatures of the medium and in turn the J/tp 
absorption in the mid and forward rapidity are approxi- 
mately the same. Therefore, the rapidity dependence of 
the regeneration becomes the key factor to control the 
experimental observation Raa^/R^a < 1 at any N p 
and (p 2 ) forward /(p 2 } mid > 1 at large N p at RHIC. 

Since the experimental results can not separate pri- 
mordial J/tp and J/tp from \ c and ip' decay, we con- 
sider a transport equation for a charmonium distri- 
bution function f^{p,x) (^ = J /tp,tp' ,Xc)- Using 
Lorentz covariant time r = \/t 2 — z 2 , space rapidity 



i] = l/21n[(t + z)/(t — z)], transverse energy Ef = 
yj E^ — p 2 with Eqi = sj m|, + p 2 and momentum ra- 
pidity ym = 1/2 In [(Ey +p z )/{Em — p z )], the three di- 
mensional transport equation can be expressed as 
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where v* = p t /Ef is the transverse velocity which leads 
to the leakage effect and is proven to be important for 
those high velocity charmonia at SPS @, and the sup- 
pression and regeneration in hot medium are described 
by the loss term a* and gain term /3$. Considering 
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the gluon dissociation process ^ + g — * c + c and char- 
monium regeneration process c + c — ► & + g in QGP, 
a (p t , y*, x t , 77, r) and /3 (p t , y*, x t , 77, r) are, respectively, 
the momentum integration of the dissociation probabil- 
ity Wg% multiplied by gluon thermal distribution f g and 

the regeneration probability W^ g * multiplied by charm 
quark distribution f c [l5| ■ is calculated with pertur- 

bative Coulomb potential 24 1 and Wff can be obtained 
from Wg^ via detailed balance. Since the experimen- 
tally measured D-meson flow is comparable with light 
hadron flow [25| . we assume charm quark thermaliza- 
tion and choose the density f c {q,x) = rA(xt)Ts(xt — 
b) cosh rj/Tda c ^ N /drjf q {q, x), where f q is the normalized 
charm quark thermal distribution, da c ^ N /drj the rapid- 
ity distribution of charm quark pairs produced in NN 
collisions, b the impact parameter, and Ta and Tg 
are the thickness functions of the two colliding nuclei 



defined as T(x t ) = lim 2 



, T(x t , zi, z 2 ) with 



T(x f ,zi,Z2) = f z 2 dzp(xt, z) and the Woods-Saxon nu- 
clear density profile /o(x). 

Since the hadronic phase appears later in the evolution 
of the fireball when the density of the system is lower 
compared to the early hot and dense period, we have ne- 
glected the charmonium production and suppression in 
hadron gas. On the other hand, from the lattice QCD 
simulation [26] on charmonium spectral function at finite 
temperature, J/tp collapses at a dissociation temperature 
T 1 / , which corresponds to the idea of sequential sup- 
pression (2?]|. Taking into account these two aspects, the 
dissociation and regeneration processes happen only in 
the temperature region T c < T < Tf of QGP, where T c 
is the critical temperature of deconfinemcnt phase tran- 
sition. Outside this region there arc a* = (3-q, = for 
T < T c and a* = /3* = 00 for T > T*. 

The transport equation has been solved analytically 
and the result is shown as 



/* (Pt,y*,x t ,?7,T) = h (Pt,y*,X*(T ),ff*(T ),T )e~ 
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with 

X*(t') = x t - vf [rcosh(j/* - rj) - t'A(t')] , 
H^(t') = y<p — arcsinh (r/r'sinh(y* — rj)) , 

A(r') = yfl + (r/r') 2 sinh 2 (y*-7 ? ). (3) 

The first and second terms on the right-hand side of the 
solution @ indicate the contributions from the initial 
production and continuous regeneration, respectively, 
and both suffer anomalous suppression in the medium. 
Since the collision time for NN interactions at RHIC en- 
ergy is about 0.1 fm/c and less than the starting time 
to of the medium evolution which is about 0.5 fm/c, the 
nuclear absorption and Cronin effect [28| for the initially 
produced charmonia have ceased before the QGP evolu- 
tion and can be reflected in the initial distribution of the 
transport equation |T]). From the PHENIX d+Au data, 
the cold nuclear matter effect is smaller at RHIC than 
those observed at lower energies |2{|. Considering fur- 
ther a finite formation time of charmonia which is about 
0.5 fm/c and larger than the collision time, the nuclear 
absorption can approximately be neglected for the calcu- 
lations at high energies [l6j]. In this case, the initial dis- 
tribution fq, (p t , y^, Xt, 77, To) is calculated according to 
the Glauber model with the J / tp production cross section 
d&% N /(ptdptdyq,) in NN collisions and the modification 
of the Cronin effect [28| . The coordinate shifts x t — > X* 
and 77 — > in the solution ([2]) reflect the leakage effect 
in the transverse and longitudinal directions. 

The local temperature, baryon chemical potential and 
fluid velocity are determined by hydrodynamic equa- 



tions and control the suppression and regeneration re- 
gion through the gluon and quark thermal distributions. 
Similar to the three dimensional transport equation for 
^> motion, a self-consistent treatment for the calculation 
of J/ip rapidity dependence needs three dimensional hy- 
drodynamic equations. Considering the fact that the ex- 
perimentally observed forward rapidity 1.2 < y < 2.2 is 
still located at the central plateau of the Bjorken hydro- 
dynamics at RHIC energy [23j , we can, for simplification, 
neglect the rapidity dependence of the fluid and take the 
transverse hydrodynamic equations [IH [l(| at y = for 
the QGP evolution at both mid and forward rapidity. To 
close the equations, we take the equation of state [3(| 
of ideal gases of partons and hadrons with a first order 
phase transition at T c , and the initial condition for the 

hydrodynamics at RHIC is the same as in Ref dam. 

We now fix the charmonium and charm quark distri- 
butions in NN collisions. With the precise transverse 
momentum and rapidity distributions for J / ip measured 
in proton-proton collisions at RHIC energy [3l|, the dis- 
tribution d<7* N I {ptdptdyy) can be parameterized as 

da% N _ 2(n-l) pj ^ n da% N 

Ptdp t dyv D(y 9 ) D(y^) dyq, 

with n = 6, £>(yw) = (p 2 t ) NN (n - 2)(1 - y|/F|), 
the maximum rapidity Y# = arccosh(y / 5jv r iv"/(2m*)), 
the averaged transverse momentum square {p1)nn = 
4.14 (GeV/c) 2 and a double Gaussian distribution 
d<J%N/dyis, [31(. At LHC energy, from the CDF exper- 
iment and CEM calculation [32|, [33|], we take n = 4, 
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(p 2 t ) NN = 12 (GeV/c^ 
into account the Cronin effect 



2 fib. Taking 
jvtv is replaced by 



2 and daZ v/dm 

MAp 2 

(pf) N = (p^) NN +a(T A (x t ,~oo,z A )+T B {xt~i>,ZB,oo)). 
The second term is from the multiply scattering of the 
two gluons with nucleons before they fuse into a ^, 
and the constant a is adjusted to the data for pA col- 
lisions dni. 

For da c j^ N /dr] at RHIC, there is a large experimental 
uncertainty and the difference among theoretical model 
estimations is also significant (35|. We take a Gaus- 
sian distribution da^ N /di] = da 1 ^ N / 'dr?| j) e _T ' with 
do c ^j N I ' dr\\ „ = 120 fib which agrees with the experimen- 
tal data [3f| and a parameter 770 determined by the as- 
sumption {da% N /dr]\ v=1 . 7 ) / (da% N /dr]\ v=a ) = 1/3 which 
is in between the maximum and minimum model calcu- 
lations (35|. At LHC, we take da c ^ N /drf directly from the 
PYTHIA simulation [H. 

the 
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fractions from direct production and Xc and ip' decay as 
6 : 3 : 1, we integrated the J/ip distribution fj/^ on the 
hadronization hyper surface and calculated the J/ tp Raa 
and (pf) at mid and forward rapidity in heavy ion colli- 
sions at RHIC and LHC energies. In the upper panel of 
FigQ]we show the ratio of forward to mid rapidity Raa as 
a function of N p at RHIC. The new experimental found- 
ing of R^^/R'Xa < 1 is difficult to be understood in 
models with only initial production mechanism, since the 
suppression in forward region should be smaller than or 
at most the same as that in central region. In our calcu- 
lation with the assumption of the same medium for the 
two rapidity regions, the ratio with only initial produc- 
tion is almost a constant, the small deviation is from the 
rapidity dependence of the J /ip transverse momentum 
distribution in NN collisions, see D(y^) in (j4|). For the 
case with only regeneration, the ratio becomes a strict 
constant and is determined by the cross sections and 
CT c g in NN collisions. Since the regeneration contribution 
at forward rapidity is smaller than that at midrapidity, 
the total ratio is less than the value with only initial pro- 
duction in the whole N p region. However, it can not reach 
the regeneration limit, because at RHIC the maximum 
fraction of regeneration is only around 50% [TBI ]. 

The ratio of forward to mid rapidity (pf) is shown in 
the lower panel of FigfTJ Again it is a constant in the 
case with only initial production or regeneration. While 
the (pf) is dominated by the initial production at for- 
ward rapidity, both initial production and regeneration 
affect the midrapidity region. In addition, the initially 
produced J/ips are harder in momentum distribution. 
Therefore, the total ratio is always larger than the limit- 
ing value set by the initial production. It starts with the 
value of (p 2 )^jv ard /(p 2 }jvjv and goes up with increasing 
N p monotonously. 

While the J/ip production at RHIC is governed by the 
competition between the initial production and regenera- 
tion, the situation at LHC is very different. At LHC, the 




FIG. 1: The ratio of forward to mid rapidity Raa and the 
ratio of forward to mid rapidity (pt) at RHIC. The lines are 
theoretical calculations and the data are from [13, Gi| . 



1.5 



IS 



1.0 



IS 



1.4 



« 1-1.2 

Q. 



S 1.0 



0.8 




IN _ 

a. 



0.6 



gforward ^gmid 



Initial 



Regeneration 



Total 




100 



200 



300 



400 



FIG. 2: The ratio of forward to mid rapidity Raa , the ratio 
of forward to mid rapidity (pi) and the ratio of forward to 
mid rapidity Saa in Pb+Pb collisions at LHC. 



4 



fireball is hotter, larger, and longer lived, almost all the 
initially produced J/tpa are eaten up by the QGP. On the 
other hand, there are about ten times more charm quarks 
at LHC than that at RHIC. As a result, the regeneration 
becomes much more important and dominates the J/ip 
production in heavy ion collisions at LHC. We show our 
prediction for the two ratios in Fig(2] The midrapidity 
and forward rapidity regions are defined as \y\ < 0.9 and 
y G [2.5,4]. With increasing N p , both ratios start from 
the initial production limit and reaches the regeneration 
limit fast. At large N p , the two ratios are almost fully 
controlled by the regeneration. The rapid change of the 
two ratios at low N p is due to the important contribution 
from the initial production in peripheral collisions, like 
the case in central collisions at RHIC. 

When the J ftp regeneration becomes dominant, we 
can use the quantity Saa = Nj/^/Np to describe the 
medium effect, where Nj/^ and Njj are, respectively, 
the J r /tp number and D-meson number per unit rapid- 
ity. Since the number of initially produced J ftps is pro- 
portional to the number of binary collisions N c and the 
number of regenerated J ftps to TV,?, we have Saa ~ l/iV c 
in the initial production dominant region and Saa ~ 1 
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